Ghosts in the self-accelerating universe 



Kazuya Koyama* 

Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth POl 2EG, UK 

The self-accelerating universe realizes the accelerated expansion of the universe at late times 
by large-distance modification of general relativity without a cosmological constant. The Dvali- 
Gabadadze-Porrati (DGP) braneworld model provides an explicit example of the self-accelerating 
universe. Recently, the DGP model becomes very popular to study the observational consequences 
of the modified gravity models as an alternative to dark energy models in GR. However, it has been 
- - - shown that the self-accelerating universe in the DGP model contains a ghost at the linearized level. 

The ghost carries negative energy densities and it leads to the instability of the spacetime. In this 
^rrj , article, we review the origin of the ghost in the self-accelerating universe and explore the physical 

■ implication of the existence of the ghost. 
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INTRODUCTION 



The acceleration of the late-time universe, as implied by observations of Supernovae redshifts [1, 2], cosmic 
microwave background anisotropies [3] and the large-scale structure [4] , poses one of the deepest theoretical 
problems facing cosmology. Within the framework of general relativity (GR), the acceleration must originate 
from a dark energy field with effectively negative pressure, such as vacuum energy or a slow-rolling scalar 
field ("quintessence"). So far, none of the available models has a natural explanation. For example, in the 
(~| \ simplest option of vacuum energy, leading to the "standard" LCDM model, the incredibly small, 

^ '. PA,obs = ^ HImI < PA^theory ' (1) 

Q>^ ' and incredibly fine-tuned, 

ff^ . ^ ^ni I today I (2) 

■ value of the cosmological constant cannot be explained by current particle physics. 

On I An alternative to dark energy plus GR is provided by models where the acceleration is due to modifications 

■ of gravity on very large scales, r > H^^ (see [5, 6] and references therein). One of the simplest covariant 
models is based on the Dvali-Gabadadze-Porrati (DGP) brane- world model [7], in which gravity leaks off 
the 4D Minkowski brane into the 5D "bulk" Minkowski spacetime at large scales. The 5D action describing 
the DGP model is given by 



^ " i / ^'^^^^ + 2^ / ^^^^ ^'^^ ~ / d^^V^Cm. (3) 

The important ingredient of the model is the induced Einstein-Hilbert term on the brane (see [8] for an early 
attempt). The existence of the brane imposes the junction condition for the metric at the position of the 
brane: 

^2 

Kf^i, - Krjf^^ = — ^(7)<i^ - «4 ^G^^): (4) 

where K^^, is the extrinsic curvature and we assume the refection {Z2) symmetry across the brane. Due to 
the induced Einstein-Hilbert action, the 4D Einstein tensor appears in the junction condition. 

On small scales, gravity is effectively bound to the brane and 4D Newtonian dynamics is recovered to a 
good approximation. The transition from 4- to 5D behaviour is governed by a crossover scale Tc [7] 
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The weak-field gravitational potential behaves as [7] 



^ for r < Tc 



^ for r > r, 



The DGP model was generalized by Deffayet to a Friedman-Robertson- Walker brane in a Minkowski bulk 
[9, 10]; the gravity leakage at late times initiates acceleration - not due to any negative pressure field, but 
due to the weakening of gravity on the brane. The energy conservation equation remains the same as in GR, 
but the Friedman equation is modified: 

p + 3iJ(p + p)=0, (7) 
H^-1='-^P. (8) 

The modified Friedman equation (8) shows that at late times in a CDM universe, with p oc — > 0, we 
have 

H-^H^^-. (9) 

Since Hq > H^o, in order to achieve acceleration at late times, wc require Vc > Hq^, and this is confirmed 
by fitting SN observations [11]. Like the LCDM model, the DGP model is simple, with a single parameter 
Tc to control the late-time acceleration although the DGP model does not provide a natural solution to the 
late-acceleration problem; similarly to the LCDM model, where A must be fine-tuned, the DGP parameter 
must be fine-tuned to match observation. 

The most interesting aspect of the DGP model is that there is a possibility to distinguish the model from 
dark energy models in GR. This is because the recovery of the 4D GR on small scales is very subtle. Although 
the weak-field gravitational potential behaves as 4D on scales smaller than Vc, the linearized gravity is not 
described by GR. This is because there is no normalized zero-mode in this model and 4D gravity is recovered 
as a resonance of the massive KK gravitons. The massive graviton contains 5 degrees of freedom compared 
with 2 degrees of freedom in a massless graviton. One of them is a helicity-0 polarization. Due to this scalar 
degree of freedom, linearized gravity is described by Brans-Dicke (BD) gravity with vanishing BD parameter 
in the case of Minkowski spacetime. Thus this model would be excluded by solar system experiments. 
However, the non-linear interactions of the scalar mode becomes important on larger scales than expected 
[12-19]. Let us consider a static source with mass M. Gravity becomes non-linear near the Schwarzshild 
radius Vg = 2GM. However, the scalar mode becomes non-linear at r* = (rgr^)^/"^ (the Vainstein radius) 
which is much larger than Vg if Tc ~ Hq^. In fact, for the Sun r* is much larger than the size of the 
solar system. A remarkable finding is at once the scalar mode becomes non-linear, GR is recovered. This 
non-linear shielding of the scalar mode is crucial to escape from the tight solar system constraints. Fig. 1 
summarizes the behaviour of gravity in the DGP model (see [20] for a review on the DGP model). 

The fact that the linearized gravity is not described by GR offers an exciting possibility to distinguish this 
model from the dark energy models in GR. In cosmology, the modification of linear perturbations affects the 
growth rate of the structure in the Universe. Thus structure formation in the DGP model is different from 
that in GR even if the expansion history of the Universe is exactly the same. Then combining the various 
observations, it is possible to distinguish the DGP model from the dark energy models in GR. Due to this 
possibility, the DGP model has been a very popular model for modified gravity alternative to dark energy 
(for a review, see Ref. [6]). 

However, these interesting features of gravity could be signatures of a pathology of the model. The non- 
linear interaction of the scalar mode becomes important on very large scales compared with usual 4D gravity. 
If we consider a Planck mass particle, the Vainstein scale becomes = {ipir'^Y^^ ~ 1000km where ipi is 
a Planck length and rc ^ Hq^. This implies that quantum gravity corrections cannot be neglected below 
1000km [21-23]. Even if we focus on the linearized behaviour of the scalar mode, there appears a problem 
of a ghost instability [22-28]. In the case where the brane is described by de Sitter spacetime, it has been 
proved that the sc;alar mode becomes a ghost. 

In this article, we focus on the problem of the ghost that appears in the self-accelerating universe. In 
section II, we study the spectrum of linearized perturbations about a de Sitter brane. Then we identify the 
origin of the ghost in section III. For the positive tension brane, the ghost is originated from a massive spin-2 
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FIG. 1: Summary of the behaviour of gravity in the DGP model. At large scales r > Fc, the theory is 5D. On 
small scales r < rc, gravity becomes 4D but the linearized theory is described by a Brans-Dicke theory. This affects 
the laxge scale structure (LSS) and the Integrated Sachs- Wolfe (ISW) effect and its cross-correlation to LSS. Below 
the Vainstein radius r < r», the theory approaches OR. This transition can be probed by weaJi lensing and cluster 
abundance as the non-linear dynamics is important for these measures. The solar system tests also provide constraints 
on the model in the 4D Einstein phase. Prom [6]. 



graviton with mass < < 2iJ^. For the negative tension brane, the spin-0 perturbation associated from 
a fluctuation of the brane becomes a ghost. For a self-accelerating brane without tension, the ghost appears 
from the mixing between the spin-0 and spin-2 perturbations. In section IV, the spectrum with matter 
source on the brane is studied. We highlight a difference between the DGP model and the massive gravity 
model that makes it difBcult to remove the ghost by a simple modification of the model such as two-brane 
model. This ghost can be identified with the brane bending mode on small scales which is a mix of helicity-0 
components of massive spin-2 perturbations and the spin-0 perturbation. In section V, the effective action 
for the brane bending mode is discussed. The leading order non-linear interaction is identified. It defines the 
Vainstein length below which the linearized analysis cannot be trusted. We confirm again the existence of 
ghost at the linearized level. In section VI, fully non-linear solutions are discussed. These solutions indicate 
that the self-accelerating universe may suffer from instabilities even at non-perturbative level. Section VII 
is devoted to conclusions and discussions. 



II. PERTURBATIONS ABOUT A DE SITTER BRANE 



A. Background solution 

Let us consider a situation where the brane is de Sitter spacetime. The bulk spacetime is a 5D Minkowski 
spacetime and the metric is given by 

ds^ = dy^ + N{yfj^^dx''dx'', N{y) = l±H\y\, (10) 

where j^,, is the metric for the de Sitter spacetime and the brane is located a,t y — 0. The Z2 symmetry 
across the brane is imposed. The junction condition at the brane gives the modified Friedmann equation 

where a is the tension of the brane. There are two branches of bulk solutions. The solution with — sign 
is called the normal branch whereas the solution with -|- sign is called the self-accelerating solution. These 
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two different solutions correspond to different embeddings of the brane. The 4D de Sitter spacetime is a 
hyperboloid in a 5D Minkowski spacetime. In the self-accelerating universe, we take the outside of the 
hyperboloid as the bulk spacetime. On the other hand, in the normal branch, we take the inside of the 
hyperboloid. 




FIG. 2: Embedding of a de Sitter brane in a flat 5D bulk. The brane world volume is the hyperboloid in the Minkowski 
bulk. The normal branch corresponds to keeping the interior of the hyperboloid, and its mirror image around the 
brane. In contrast, for the self-accelerating branch, we keep the exterior, and its reflection. From [26]. 



B. Spectrum of perturbations 



We study the linear perturbations 



(12) 



about the background de Sitter spacetime. In addition to the gravitational perturbations /i^^, we must take 
into account a perturbation of the position of the brane y = ^{x) [30]. Using the transverse-traceless gauge 
V^/ipi, = = 0, the perturbed junction condition is given by [24] 



(13) 



where H = N'/N\y^f) = ±H, is a covariant derivative with respect to 7^,^, fc^i, = {l/2)dyh^i, on the 
brane is perturbations of the extrinsic curvature ii'^j^ and X^i, is given by 



(14) 
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The equation of motion for ip is obtained from the traceless condition /i = 0; 

(l-2Wre)(a4 + 4if>=^, (15) 

where ^4 = V^V^. 

Let us find solutions for the vacuum brane T^i^ = 0. Using the separation of variables h^v = 
/ dm env{x)um{y), the equation of motion in the bulk is written as 

< + ^(m2-2if>„=0, [a^ - - 2H^)e^,ix) = 0, (16) 

where prime denotes a derivative with respect to y. There are two types of solutions to Eq. (13). One type 
of the solution is a homogeneous solution with <fi = 0, which is called the spin-2 perturbation. The spin-2 
perturbations Xju^ satisfy the junction condition without 

'X'nu - '^'X-Xnu = -m\cXnu- (17) 

We find a tower of continuous Kaluza-Klein (KK) modes starting from = (9/4) if ^ as well as a normal- 
izable discrete mode. In the self-accelerating branch, the solution for the discrete mode is given by 

n„,«iV(y)-i+i/«'-=, g| = ^^(3i?re-l), (18) 

for Hrc > 2/3 [29]. For Hrc > 1, the mass is in the range < < 2H'^ where = 2H'^ for the 
self-accelerating universe Hrc = 1 and — > for Hrc — * 00. In the normal branch, the discrete mode is a 

zero mode = 0. 

The other type of solution to Eq. (13) is an inhomogoncous solution sourccd by the scalar mode (p. We 
call this solution the spin-0 perturbation. In the self-accelerating branch, there is a normalizable solution 
given by [24] 

^ H{l-^Hrc) ^^"^^ ^ (^^^ 

This is a solution with = 2H^. In the normal branch, there is no normalisable spin-0 perturbation. Then, 
the normal branch is ghost-free. Eq. (19) is singular at Hrc = 1 and we will deal with this case separately 
in section III.C. 



III. GHOST IN DE SITTER SPACETIME 
A. Effective action 

We can construct the 2nd order action for h^i, and (fi from the 5D action (3). The result is given by [24] 

S2S = J d^x^f^N-^h''''5^^^G^. + d^x^f^CB, (20) 

where d'^^^G^^ is the 5D perturbed Einstein tensor and 

Cb = k^"'h^,-kh+^nih'' -h^^h^,) 

+ {1- 2Hrc) {h^uVV'fi - hVfVpip - SH'^hcp) 

- 3n(^-{l-2nrcMa^+iH^)ip + YT^^ 

+ ^K^h'^-T^^-!lh^^-X^4h). (21) 
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This action gives the correct equation of motion and the junction condition for h^i, and the equation of 
motion for (p. 

In the following discussions, wc focus on the sclf-accclcrating branch. We can derive an effective action 
for the brane fluctuation (p by substituting the 5D solution for h^^ given by ip (19) into the 5D action (20) 
and get the off-shell action for (p by integrating out only with respect to the extra coordinate y [31]. This 
yields the action for ip as 

The 4D effective action for the spin-2 perturbations is also obtained in a similar way. For the discrete mode 
with m^, we get 

' 4%^; - 2) / '^"^X-iu. - - (23) 

where transverse-traceless gauge fixing conditions V^x^i^ — X'^ — ^ Eire imposed . This is exactly the same 
action for the spin-2 perturbations in the 4D massive gravity theory where the Pauli-Fierz (PF) mass term 
is added to the Einstein-Hilbert action by hand [32] 

Sm = / d^xV^ih^'-'h^,^ - h^). (24) 



B. Ghost in de Sitter spacetime with a tension 

Ref. [24] studied the existence of the ghost based on the above effective action. In the limit Hvc — > oo, the 
action (23) approaches the one for massless spin-2 perturbations. However, there is a discontinuity between 
the massless perturbations and the massive perturbations that is known as the van Dam-Veltman-Zakharov 
discontinuity [33]. Due to the lack of gauge symmetry, the massive spin-2 perturbations contain a helicity-0 
excitation. Moreover, it has been shown that this helicity-0 excitation becomes a ghost if < < 2H^ 
[34, 35] (see also [26]). This is exactly the same mass range for the discrete mode in the self-accelerating 
branch for Hvc > 1. Thus we identify the ghost as the helicity-0 mode of the discrete mode of the spin-2 
perturbations. On the other hand, for Hrc < 1, the spin-2 perturbations become healthy as the mass for the 
spin-2 perturbations is larger than 2H^. However, in this case, the coefficient in front of the effective action 
for the spin-0 perturbations becomes negative and the spin-0 perturbation becomes a ghost. 
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FIG. 3: Summary of the mass spectrum of the scalar perturbations in + branch. Spin-2 perturbation has continuous 
modes with > (9/4)// and a discrete mode with rn^ = while spin-0 perturbation has = 2H^. In the limit 
Hrc — » 1, both the helicity-0 excitation of spin-2 perturbation and the spin-0 perturbation have mass m = 2H^ and 
there is a resonance. 
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FIG. 4: Summary of the existence of the ghost. From [24]. 



C. Ghost in self-accelerating universe 



In the sclf-accclcrating universe Hrc = 1 , the mass of the discrete mode of the spin-2 perturbations becomes 
2H^. This is a special mass in the massive gravity as the action is invariant under the transformation 

+ (V^V, - H^-f^,)X, (25) 

where X is any solution of the equation (^4 + AH'^)X = 0. This is so called enhanced symmetry and the 
helicity-0 mode can be eliminated by this symmetry [35]. Then there are only four polaization when the 
mass is = 2H^ and the ghost disappears. 

However, this does not happen in the self-accelerating universe as there is an additional spin-0 perturba- 
tions with the same mass [25]. In fact in the self-accelerating universe, Hrc = 1, the bulk wave function for 
the discrete mode is given by u„i,f 3c N{y) and it is the same as the spin-0 perturbation. Then the discrete 
mode of the spin-2 perturbations and the spin-0 perturbation degenerate and they can mix. This comes 
from the fact that at = 2H^, the spin-0 and spin-2 perturbations degenerate. In fact we can make the 
spin-2 perturbations from a scalar 

h^^f^ = {V^V.-H^j^,)X, {n^+4H^)X = 0. (26) 

This is a scalar mode with mass squared —AH^. At the same time /ij^i/ is a transverse-traceless perturbations 
and from the identity, 

(□4 - 4if2)/j(^2H^) = (v^V. - if V) (°4 + ^H^)X = 0, 

it is identified to have a mass given by = 2H^. 

Let us investigate the limit Hrc — * 1 carefully. The solutions for the metric perturbations are given by 



1 — 2Hr 



(27) 



where is given by Eq. (18) and we only consider the localized modes. The limit Hrc 1 looks singular, 
but we can perform a field re-definition 

X™<'(a;) = A,, - -Azi^(v^V. + H^7,.)^. (28) 

The metric perturbations read 

1 — 2Hr 

h^,{x,y) = A^^UraM + _ Hrc) ^^"^^ ^ ^^^f^-M^ ' ^mAy))- (29) 
Then taking the limit Hrc — > 1, we get [25] 

V = ^f^-i^) + ;^(V,.V. + H^j^^)^{x) log (1 + Hy) . (30) 
Substituting this expressions into the bulk equation and the junction condition, we get an equation for 

□4^^. - 4H^A^^ = H (V^V. + >fi. (31) 
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The effective action for A^^, and (p is obtained by substituting the solutions for A^i^ and (p into the 5D 
action and integrating out the extra-dimension as 

Seff = ^J d^xV^^ - A^'-X^^A) - H^A^'^A^, + H^A^ 

where we introduced the notation A = A'^. The first hne in Eq. (32) coincides with the quadratic Lagrangian 
for the Pauli-Fierz theory of massive gravity with = 2if^. Thus the action could have the enhanced 
symmetry (25) if there were no mixing to However, the mixing between and A^^fj, breaks this symmetry 
explicitly. 

In order to study the existence of the ghost, it is necessary to write the effective action only in terms of 
physical degrees of freedom. Ref. [25] performed the Hamiltonian analysis to derive the reduced Hamiltonian 
written only in terms of physical degrees of freedom. There are two dynamical degrees of freedom, namely 
the spin-0 mode and the helicity-0 excitation of spin-2 perturbations. It is found that, in general, the 
Hamiltonian cannot be diagonalized. However, on small scales under horizon, it is possible to diagonalize 
the Hamiltonian and we find a ghost from a mixing between the spin-0 perturbation and the helicity-0 
excitation of spin-2 perturbations. 

D. Two-brane model 

Ref. [28] tried to remove the ghost by adding the second brane in the bulk. In previous section, we saw 
that the origin of the ghost for the positive tension brane is coming from the fact the discrete mode for 
the spin-2 perturbations has a mass in the range < m? < 2H^. It is possible to eliminate this ghost by 
introducing the second brane in the bulk. By making the distance between the two branes short, the mass 
increases as in the usual Kalzua-Klein compactification. Once the mass becomes larger than 2H^, the spin-2 
perturbations do not contain a ghost. In fact for Hvc = 1, once the second brane is put and it cuts ofi^ the 
bulk spacetime, the mass becomes > 2i?^ regardless of the distance between the two branes. Note that 
even in the presence of the second brane, the Friedman equation on the visible brane is unchanged and the 
brane can self-accelerate. 

However, as soon as the mass of the discrete mode of the spin-2 perturbations exceeds = 2H^, it is 
proven that the spin-0 perturbation becomes a ghost [28] . At the critical length between the two branes where 
the mass of the discrete mode of the spin-2 perturbations is given by = 2H^, the spin-0 perturbations mix 
with the helicity-0 component of the spin-2 perturbations exactly in the same way as the self-accelerating 
universe in one brane model and there appears a ghost. 

The spin-0 perturbation is the radion that describes the distance between two branes. Then one would 
try to remove this ghost by eliminating the radion by stabilizing the distance between two branes. Once 
we stabilize the radion, the brane fluctuation mode with mass = 2H'^ becomes non-physical. The 
simplest way to achieve the stabilization is to introduce a scalar field in the bulk. Then instead of the brane 
fluctuation, there appears scalar field perturbations which can have a discrete mode and an infinite ladder 
of massive modes. In general, the mass of the discrete mode is different from — 4iJ^ which corresponds to 
= 2H^ in the spin-2 perturbations language. However, it turns out that if the mass of the discrete mode 
of the spin-2 perturbations becomes = 2H^, the discrete mode of the scalar field perturbations has a 
mass = —4H^ which is the special case where it can mix with the spin-2 perturbations. Then completely 
the same phenomena happens as in the one brane model. Once the mass of the discrete mode of spin-2 
perturbations exceeds 2H^, the helicity-0 mode of the spin-2 perturbation becomes non-ghost but the scalar 
field perturbation becomes a ghost! Then it is impossible to remove the ghost. In the next section, we will 
explain why it is impossible to remove the spin-2 ghost and spin-0 ghost simultaneously. 

IV. SPECTRUM WITH MATTER SOURCE 

In the previous section, we identify the origin of the ghost by the analysis of the spectrum without matter 
source. The most interesting finding is that it is impossible to remove the spin-2 ghost and the spin-0 ghost 
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simultaneously [26, 28]. We can clarify the reason for this by studying the spectrum with matter source. We 
also derive the effective theory for perturbations on small scales. 



A. Amplitude 



Wc introduce matter source on the brane [26-28]. The junction condition for the transverse- traceless 
modes inchiding matter perturbations is given by Eq. (13). Using the equation of motion for the brane 
bending Eq. (15), Eq. (13) is given by 

{dy - 2H) /^(^) = -K^S,. - re(D4 - 2if^)/.(r), (33) 

where 

= V - lH^7^^-T + i(V^V. + ff27M-)(°4 + 4ij2)-ir. (34) 

The solution for transvcrsc-tracclcss perturbations can be obtained by the Green's function which can be 
constructed by solutions without source obtained in section II. The solution is written as 

i * 

Here the solution for the properly normalized mode functions Ui(0)^ are given by 

2,„, 1 3Hr,- 2 2 3ffrc-l 



«m(0) = -1 Z-2 ' (36) 



We should bear in mind that in the TT gauge, the brane is not located at y = and the brane bending must 
be taken into account. The induced metric on the brane is given by 

^in<i„ced(o) ^ ^TT _ 2H^^,^. (37) 

The one particle exchange amplitude is defined as 

A = ^h'|^^'^'''"^{Q)T^"' = ^h'-^J^T^"' - HipT. (38) 



Using Eq. (35) and Eq. (15), A is calculated as 
A = -K^Y^UiiO) 



- 2H^ - mf 3 □4 + 6^f2_^2- 



~l-2Hrc~6^ai + Am^- ^^^^ 

The second line becomes a double pole when m? = 2H^ which occurs in the self-accelerating universe 
Hvc = 1 [27]. Usually, a double pole can be recast into difference of two simple poles, giving rise to a ghost. 
In fact, we can rewrite the second line as 
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The first term has a pole at 04 = — 6i?^ which corresponds to the spin-2 contributions. The second term 
has a pole at ^4 = —AH"^ which corresponds to the spin-0 perturbation. Around m? = 2iJ^, these terms 
become dominant in the amplitude and determine the existence of the ghost. A crucial point is that the 
contribution of the spin-0 perturbation is always opposite to that of the spin-2 massive perturbations [28] . 
We have already seen that the massive spin-2 perturbations contain the helicity-0 mode that is a ghost if 
< m? < 2H'^. In fact we see that the sign of the amplitude changes at m? = 2iJ^. Then for < < 2i?^, 
the spin-2 perturbation mediates the repulsive force. This means that spin-0 perturbation is not a ghost as it 
must mediate the opposite force compare to massive spin-2 perturbations. On the other hand, for 2H'^ < mf , 
the spin-2 perturbations do not carry a ghost and mediate a normal force. Then the spin-0 perturbations 
should mediate a repulsive force and it should be the ghost. This explains the reason why the spin-0 ghost 
appears as soon as the spin-2 ghost disappears. 

In the self-accelerating universe, there appears a double pole [27]. This manifests the fact that the spin- 
2 and spin-0 perturbations degenerate. In fact if we try to separate the poles for the spin-0 and spin-2 
perturbations, we encounter a divergence. Then we need a careful treatment for this special case. 

It should be emphasized that the existence of spin-0 contribution is crucial to have a non-singular ampli- 
tude. It is instructive to compare the amplitude (39) with that in massive Pauli-Fierz theory [36] 



"^massive 



1 lM2-3i?2 

Ta,.- TTTT^ TT^T^"" - - —5 — 



□ 4-2il2_M2 3M2-2iJ2 □4-^6if2_M2" 



(41) 



We can check that the spin-2 contributions in (39) are the summation of these massive spin-2 perturbations. 
The amplitude is singular for = 2H^. This comes from the fact that in Pauli-Fierz theory, we cannot 
couple gravity to matter with non- vanishing trace of energy momentum tensor. This is because the equation 
of motion is given by [36] 

{2H^ - M2)(n4 + 4H^)h = ^^^T, (42) 

where h is the trace of perturbations and T must vanish for = 2H^. This is the origin of the singularity 
in the amplitude and this is a pathology of massive gravity theory in de Sitter spacetime. In the DGP 
model, it is the brane bending mode that couples to the trace of energy momentum tensor and there is no 
such a restriction even if the discrete mode has a mass mf = 2H^. In fact the singularity at mf = 2H'^ is 
exactly canceled by the spin-0 perturbations as can be seen in Eq. (40). Thus in the DGP, we do not have 
the pathology and we can couple gravity to matter perturbations with non-zero trace of energy- momentum 
tensor. However, we have seen that it is this mechanism that makes the spin-0 perturbation a ghost when the 
spin-2 perturbation becomes non-ghost: in order to cancel the singularity in the spin-2 sector at m| = 2H^, 
the spin-0 contribution must be opposite to that of the spin-2 perturbations. 



B. Effective theory on small scales 



On small scales, it is possible to derive the 4D effective theory for perturbations. Let us consider the limit 

□4>i^^m^ (43) 

Then the amplitude is approximated as 



1 



1„ 1 



Tut. — Tf"" - -T — T 



□ 



3 04 



H 



T—T. 



6 1 - 2Hrc ^4 



(44) 



Note that the self-accelerating universe Hrc = 1 is not a special case anymore. Now using the solutions for 
the mode functions we get 

Y,MO)' = «d(0)'+^zi,„(0)2 

i 

1 3Hrc-2 H 

dk- 



2re 3Hrc - 1 n{Hrc)^ Jo {k"^ + |) (A;2 + 9 _ ml) 
1 

2^r' 



(45) 
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The effective gravitational coupling is read as 



(46) 



Thus the 4D gravity is recovered by the summation of massive gravitons. Then the amplitude is calculated 
as 



-K4 



1 



1„ 1 



Tnu — T'*'' - -T — T 



□ 



3 D4 



y 1 - 2Hrc' ^4" 



-T — T. 



(47) 



The first terms is exactly the same as the amplitude of massive gravity in the Minkowski spacetime. Due to 
the scalar polarisation, the coefficient in front of TU^^T is 1/3 not 1/2. The last term represents the effect 
of the curvature of the brane. Then finally wc get 



A 



1 



"''u/ 31- 



ll-3iJr^„ 1 



-T — T 
2Hrc 



(48) 



This result can be compared with the 4D Brans-Dicke (BD) theory. In the BD theory with BD parameter 
CO, the amplitude is given by 



A = -kI 

Then the BD parameter is determined as 



□4 



31 + ^ 04 



u = —ZHvc 



(49) 



(50) 



This agrees with the results obtained in Refs. [16, 17, 37, 38]. It is known that the BD theory contains a 
ghost if 

■i (51) 



This means that there is a ghost if Hvc > 1/2, which agrees with the spectrum analysis. 



C. Boundeiry effective action 



The analysis in the previous section indicates that it is possible to derive the 4D effective action on small 
scales. Let us again consider the limit 

□4 > r-2, H'^ (52) 

Using this approximation, wc can only keep the 4D terms in the 5D second order action (20) and (21). Then 
the 4D boundary effective action is obtained as 



(1 - 2Hrc) {h^.^VVcp - hWWp^ - SH'^h'fi) 
,^2 



-m { -(1 - 2nrc)<p{ai + AH')^ + —T^ 



(53) 



where Ti = H in the self-accelerating branch. It is possible to check that this action consistently reproduce 
the amplitude (48). We can diagonalise the action by defining 



The resultant action is 



(54) 



(55) 



The sign of the kinetic term changes at Htc = 1/2 and we find that the brane bending mode becomes 
a ghost for > 1/2. 
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V. BOUNDARY EFFECTIVE ACTION 

In the previous section, we show that it is possible to derive the 4D effective action. Refs.[22] and [23] 
derived the 4D effective action for the brane bending mode including the non-linear interactions and identified 
the scale at which the linearized analysis cannot be trusted. 

A. The effective action for the brane bending 

Refs. [22] and [23] derived the boundary effective action for perturbations on Minkowski background, 
guN = rjMN + huN- They integrated out the bulk to obtain an effective action for the 4D field living on 
the boundary. The final result at the quadratic order is 



(56) 



where A is a non-local differential operator, A = ^/—04 and the kinetic terms have been diagonalized by 
defining 

hyy = -2A'!r, n^ = Nf^-d^TT, x^i^ = - r^^TTTj^^, (57) 

and is a shift which is (y, /u)-component of the 5D metric, gy^. In this gauge, tt plays the same role as 
the brane bending (p. Note that this result is consistent with the previous effective action (55) for tensor 
and scalar parts if one takes H = 0. 

By taking into account bulk interaction with higher powers of Hmn, one finds that the leading order 
boundary interaction terms is cubic in tt, and involves four derivatives. 



5(3) 



~ J d'^xidTrfa^TT. (58) 



In order to extract the non-trivial non-linear interactions, let us consider the case where the flat approxi- 
mation is good Xftv *C 1- However we want to preserve the self-coupling of the tt field. In terms of canonically 
normalized field tt defined by ^ = 7r/(2K4rc), the cubic self-coupling is unchanged if we keep 

/2Miy/^ _ /M^y/^ j_ j_ 



fixed (see below). We also preserve the interaction between n and the energy momentum tensor The 
interaction between hfj^u and T^^ is given by {l/2)hf^,^T^'^ . Prom the definition of tt, n interacts with matter 
via {1 / M4)ttT . Therefore, if we take the formal limit 

M4 — > 00, Tc — > 00, T^i, — > 00, A = const, — ^ = const. (60) 

we can decouple 4D gravity while keeping the full Lagrangian for n. It is possible to check that all further 
interactions other than the cubic interactions vanish. Then the full action for the tt field in the flat spacetime 
is given by 



3(S7r)2--l(a7r)2D47r+^^r 



From this action, the equation of motion for tt is derived as 



(61) 



3°47r - j^id^d^nf + ^(□47r)^ = (62) 

It is possible to understand this equation from a geometric point of view. The Gauss-Codazzi equation in 
the bulk is given by 

(^)i? + Ki.^K^'' -K'^ = 0. (63) 
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The junction condition at the brane is given by 



K^,. - Kt]^. = -^{T^, - KfG^^). (64) 



Then combining these equations we get 



Note that this equation is exact as this is a combination of a geometric identity (63) and the junction 
condition on the brane (64). Prom Eq. (57), the extrinsic curvature is calculated as 

Kfj,t, = ^dij,d^n, (66) 

in the limit (60) . Then Eq. (65) gives the equation of Eq. (62) . 

Once a solution of Eq. (62) or Eq. (65) is found, we can perturb tt and expand the action up to quadratic 
order in the perturbations p. We find 

6S = j d^x [-^dpf - 2{K^^ - i^^^k)d^pd''p^ . (67) 
where K^^, = TcK^^, and it satisfies 

3^ + ^,.^--^^ = -^. (68) 

B. Self-accelerating universe 

Let us consider a de Sitter background solution with tension a. The equation for K^^ admits two branches 
of solutions: 



This corresponds to de Sitter solution where the Hubble parameter H is given by 

H = \n\ (70) 

The self-accelerating branch corresponds to + sign in Eq. (69). It is also possible to derive a corresponding 
solution for ft. The de Sitter invariant solutions are 

fro = -t^A^xf^x^. (71) 

Now let us consider a small fluctuations around the de Sitter spacetime. The action for the fluctuations 

p (67) is given by 

5S = - j d'^x [3 (1 - 2Wrc) {dpf] . (72) 

This agrees with the result for perturbations around de Sitter solution in the self-accelerating branch in the 
limit ^4 ^ H^. In particular, the fluctuation p becomes a ghost when Hvc > 1/2 in the self-accelerating 
universe. In fact this was the first discovery of the ghost in the self-accelerating universe in literatures. Note 
that for large Tire, the kinetic term becomes large and the scalar mode becomes non-dynamical. In the 
BD theory, this corresponds to the limits where the BD parameter u = —ZTirc becomes infinite. Then we 
recover GR at the linearized level [39]. In this limit, the non-linear interactions are also suppressed. 
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C. Non-linearity interactions 



Next, wc study the effect of tlic uou-Iiiiear interactions of ir. Let is consider a static localized soiu'cc of 
mass M. We look for a static spherically symmetric solution ^{r) [22, 23]. The equation for n (62) is then 
given by [18, 19] 



^ + ^^.^3^^ + ^) = ^ (73) 
dr^ rdr J ^ ^ ' 2M4 ' ^ ' 



where /? = 1 — 2Hrc and 



„ _ 2 
"~ A3 



Then it is possible to integrate the equation to get 



3/37r + - + ^^ = 0, (75) 



where 



is the Schwarzschild radius of the source. Hereafter, we assume Vg — const, for simplicity. Taking the r 
derivative of Eq. (75) gives an algebraic equation for dn/dr. Then we get a solution for dn/dr as 



where 

which is the Vainshtein radius for a source. On scales larger than r*, the non-linear interactions can be 
neglected. On the other hand, on scales smaller than r, the non-linear interactions cannot be neglected and 
the linearized analysis cannot be trusted. Note that r, is much larger than Vg if Tc ^ Hq where Hq is the 
present-day horizon scale. 

D. Spherically symmetric solution 

Let us study the behaviour of the spherically symmetric solution [14, 16-19]. The metric perturbations 
are given by 

2 

where Xixu satisfies the 4D GR equations. We take the scalar perturbations on a brane as 

ds"^ = -(1 + 2*)dt2 + e2^*(l + 2^)5ijdx'dxK (80) 
Then the solutions for the metric perturbations can be obtained as 
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On scales larger than the Vainshtein radius r > , the solutions are given by 

These solutions can be described by BD theory with BD parameter w = 3(1 — (3)/2 = —SHrc- This agrees 
with the analysis in the previous section. On scales smaller than the Vainshtein radius, r < r*, the cubic 
interaction cannot be neglected. Then the linearized analysis cannot be trusted. In this case, the solutions 
for \1/ and $ are obtained as 

$ = J+sign(/3),/:^, (85) 

* = -r^+sign(/3)J^. (86) 

In this region, the corrections to the solution in 4D GR are suppressed for r < r* so that Einstein gravity is 
recovered. Prom Eq. (73), we can see that S dominates over the linear term in this region. This indicates 
that once tt becomes non-linear, the solutions for the metric approach those in 4D GR. We should note that 
/3 is negative in the self-accelerating solution while [3 is positive in the normal branch solution. Then the 
corrections to 4D GR solutions have opposite signs in these solutions, as was first pointed out in Rcf [16]. This 
means that the metric perturbations on small scales r <C is sensitive to the cosmological background 
solutions. We should note that Ref. [40] claimed that these solutions do not satisfy the full set of the non- 
linear equations. However, Ref. [19] showed that these solution are fully consistent with the 5D non-linear 
equations as long as we consider scales larger than rg. 




VI. NON-PERTURBATIVE SOLUTIONS 



In the previous section, we saw that the non-linear interaction of the brane bending mode is important 
below Vainstein length r*. This also means that as long as we consider length scales larger than Vainstein 
radius, the linearized analysis can be trusted and we find a ghost. However, in massive gravity theory, it is 
known that the non-linear interactions are very subtle. As in the DGP model, the perturbative approach 
that takes into account non-linear interactions of the scalar mode (helicity-0 excitations) shows that the 
solution approaches GR near the source [41, 42]. However numerical attempts to find fully non-perturbative 
solutions have failed in massive gravity theory [43]. Then it is required to study fully non-linear solutions 
carefully to check the validity of the linear perturbations [44] . 



A. Schwarzschild solution 



The attempt to find fully non-perturbative Schwarzschild solution was made in Ref. [45]. They assume 
that the 5D metric takes the form 

ds^ = -e'^dt'^ + e^dr^ + r'^dQ'^ + 2'ydydr + e^dy^, (87) 

where y is the extra-dimensional coordinate. In general, we should solve 5D Einstein equations and impose 
the junction conditions to find solutions. Instead of solving the bulk, they assumed u = —A to close the 
equations without solving the bulk equations. Then, remarkably, they can manage to find analytic solutions 
for A. In the normal branch, below the Vainstein radius, the solution is given by [45, 46] 

r yrcJ^rJ 

which is similar to the solution obtained in the previous section. However, on scales larger than Vainstein 
scale, they find 

/ \ 1/3 
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which is the 5D solution and quite different from the previous 4D BD solution. Especially the mass of the 
source for observers located at r > r# is 'screened'. In a naive perturbative approach, we would expect that 
the large scale Newton potential behaves as 

-^ = -^' 5 = 2G5M = 2rerg. (90) 

In the solution (89) the 5D gravitatinal length is smaller fg <C r^ ^. This indicates that the effective mass 
M measured at r > r* is screened by M = M{rg/rcY/^. 

If this solution is a true non-perturbative solution in the DGP model, this means that the linearized 
perturbations cannot be trusted even at r > . However, we should bear in mind that the solution is 
obtained by closing the equations by an ad hoc metric ansatz. It is far from trivial that the solution obtained 
in this way satisfies the reasonable boundary conditions in the bulk. In fact it is pointed out that the 
regularity condition in the bulk is crucial to specify the behaviour of gravity on the brane [38] . It is an open 
question to check that the solution obtained in Rcf. [45] has a regular bulk solution. 

However, this solution also shows interesting features in the self-accelerating background that might be 
related to the ghost. In the self-accelerating universe, the effective mass measured at r > r* becomes 
negative. This suggests that the self-accelerating background may not be problem-free even in the full 
non-linear theory. 



B. Domain wall 



The other example of the exact non-linear solutions is a domain wall. Ref. [47] found an exact domain wall 
solution. In the 4D spacctimc, the domain wall creates a jump in the extrinsic curvature. On the other hand, 
the domain wall in the 5D spacetime is a co-dimension 2 object like a string. Then the domain wall creates 
a deficit angle in 5D spacetime. The relation between the deficit angle 5 and the tension of the domain wall 
fjb is given by 

2ffretan7±7=^, (91) 

where 7 is related to the deficit angle as (5 = t47, where — is for the self-accelerating branch and -|- is for 
the normal branch. The first term comes from the jump of the 3D extrinsic curvature in the brane as in 4D 
GR and the second term is the contribution from the deficit angle in the 5D spacetime. 
For small tension a <C M5 , the deficit angle is given by 

5= ^- ^, n = ±H. (92) 

In the normal branch, there is a screening of the tension. The deficit angle in the 5D spacetime is smaller 
than expected. On the other hand, in the self-accelerating branch solution, there is 'over-screening'. This 
indicates that the wall in the self-accelerating universe behaves as with negative tension. Then it appears 
that the stability of the solution is not granted. This may be related to the existence of the ghost. In fact, we 
see that the factor that determines the screening 1 — 2'Hrc is exactly the factor that determines the existence 
of the ghost. 



VII. CONCLUSIONS AND DISCUSSIONS 



In this review, we discuss the ghost problem in the self-accelerating universe. First we studied the spectrum 
of the perturbations without matter perturbations about a de Sitter brane. For the positive tension brane, the 
ghost comes from the fact that there is a discrete mode for the spin-2 perturbations with mass < < 2iJ^. 
The massive spin-2 perturbations contain a helicity-0 mode that becomes a ghost if the mass is in the range 
< < 2H^. In the self-accelerating universe without tension, the mass becomes = 2H^. This 
is a special mass in Pauli-Fierz massive gravity theory because there exists an 'enhanced symmetry' that 
eliminates the helicity-0 mode. However, in the DGP model, there is a spin-0 perturbation with the same 
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mass and this breaks the symmetry, leading to a ghost from the mixing between the spin-0 and spin-2 
perturbations. For a negative tension brane, the spin-0 perturbation becomes a ghost if Hr^ > 1/2. It is 
easy to remove the spin-2 ghost by putting a second brane in the bulk and make the distance between the 
two brancs small. Then the mass of the discrete mode of the spin-2 perturbations becomes larger than 2H^ . 
However wc find that the spin-0 perturbation, the radion, becomes a ghost. If we stabilize the radion, the 
perturbations of the scalar field that is necessary to stabilize the radion becomes a ghost. 

We recovered the same result by studying the spectrum with matter source. The one particle exchange 
amplitude is the summation of massive spin-2 perturbations and a spin-0 perturbations. The amplitude of the 
massive spin-2 perturbations diverges if the mass is = 2iJ^. In the DGP model, the spin-0 perturbation 
exactly cancels this singularity at Htc = 1 where the spin-2 mass becomes m? = 2H'^. This means that 
the spin-0 interaction is opposite to that of massive spin-2 perturbations. This is the reason why the spin-0 
perturbation becomes a ghost once the massive spin-2 perturbation becomes healthy. At small scales, the 
effective theory for perturbations is described by the BD theory with BD parameter w = —3Hrc- Here the 
BD scalar is a mix of the spin-0 perturbation and the helicity-0 component of spin-2 perturbations. In the 
BD theory, the BD scalar mode becomes a ghost if w < —3/2 and this condition is given by Hrc > 1/2 which 
agrees with the spectrum analysis. This scalar can be identified as the brane bending mode. It is possible 
to construct the effective action for the brane bending by keeping the boundary terms in the full 5D action. 

The effective action for the brane bending was shown to be a powerful tool to analyze the non-linear 
interactions of the scalar mode. In the decoupling limits where gravity decouples from the scalar interactions, 
the cubic interaction is the dominant non-linear interaction. For a local source with the gravitational length 
Vg, the scalar mode becomes strongly coupled below the Vainstein length r* = (r^r^)^/^. Then the linearized 
perturbations cannot be trusted below this length scale. Wc checked that the effective action for the brane 
bending reproduces the results obtained in the linearized analysis beyond r*. Especially, the brane bending 
mode becomes a ghost in this linear regime if Hrc > 1/2. 

This seems to indicate that the linearized analysis makes sense at least at large scales r > and we 
cannot avoid the ghost. However, in massive gravity theory, fully non-pcrturbative effects are very subtle 
and it is possible that the linearized solution cannot be matched to the fully non-perturbed solutions. There 
are two known fully non-perturbed solutions. One is Schwarzschild solution obtained in Ref. [45]. Although 
it is an open question that this solution has a physical solution in the bulk, the solution does show that the 
effective mass measured at r > becomes negative in the self-accelerating universe. This may be related to 
the ghost in the perturbative solution. The other solution is a domain wall solution. In this case, the exact 
5D solution is known. For a small domain wall tension, the domain wall has a negative deficit angle from 5D 
point of view in the self- accelerating branch solutions. This is yet another evidence that the self-accelerating 
universe suffers from the ghost instability even at non-perturbative level. 

There are many open questions that deserve further studies. We will address some of issues. 

Instabilities 

The ghost has a wrong sign for its kinetic term. This means that the energy density can be indefinitely 
negative. This can lead the classical instability of the system. In Ref. [26] it is argued that the self-accelerating 
universe must be classically unstable. However, the ghost is found at linearized level and so far no classical 
instability was found at linearized level. We should carefully study the non-linear interactions to address 
the classical stability of the model and this is still an open question. See Ref. [48] for recent discussions on 
non- linear instabilities. 

Quantum mechanically, the ghost leads to spontaneous pair creation of ghosts and normal particles. Once 
such a channel opens, Lorentz invariance leads to a divergence of the particle creation rate and the decay 
rate of the vacuum is infinite. The same problem occurs in so-called phantom cosmology where the ghost 
scalar field is introduced to explain the dark energy with equation of state smaller than —1. It is argued 
that we can avoid the rapid decay of the vacuum if there is a Lorentz non-invariant UV cut off of the order 
MeV in the phantom cosmology [49]. In the DGP model, the situation is more subtle [50]. If we consider 
the situation in which there is a spin-2 ghost, we need to treat the helicity zero mode in a different way. 
Otherwise, negative norm states appear. But if we take a different prescription for the quantization for the 
helicity zero mode, this procedure necessarily breaks de Sitter invariance. When there is a spin-0 ghost, a 
similar phenomenon happens. In this case, the mass of the ghost is given by —4H^. But we know that 
there is no de Sitter invariant vacuum state for a scalar field with negative mass squared. Once de Sitter 
invariance is broken, one may be allowed to consider the possibility that the non-covariant cutoff scale may 



18 



arise due to the strong coupling effect. Tfie strong coupling length is very large ~ 1000 km. Then the 
particle creation could be milder than the usual ghost in the Minkowski background [51]. Clearly, further 
studies are necessary to verify this. 

The fate of instabilities 

If the self-accelerating branch solutions have the ghost instability, then we are naturally lead to ask what 
does this solution decay to. An interesting fact is that the normal branch solutions are ghost-free. Then 
it is tempting to think that the self-accelerating solutions decay into the normal branch solutions. In fact 
for a given tension, the Hubble parameter in the self-accelerating universe is larger than that of the normal 
branch solutions. Then there could be the nuclcation of bubbles of the normal branch in the environment of 
the self-accelerating branch solution. This would resemble a kind of false vacuum decay in de Sitter space. 
False vacuum decay is described by an instanton which is a classical solution in an Euclidean time connecting 
initial and final configurations. In our case, we are interested in a solution that interpolates between the self- 
accelerating and the normal branches. Ref. [50] tried to construct such an instanton solution. It was found 
that the solution requires the presence of a 2-branc (the bubble wall) which induces the transition. However, 
this instanton cannot be realized as the thin wall limit of any smooth solution. Once the bubble thickness 
is resolved, the equations of motion do not allow 0(4) symmetric solutions joining the two branches. It was 
concluded that the thin wall instanton is unphysical, and that one cannot have processes connecting the two 
branches. This suggests that the self-accelerating branch does not decay into the normal branch by forming 
normal branch bubbles. Thus it is still unclear what is the end state of the ghost instability. 

Possible ways out of the ghost? 

In Ref. [27], several ways out of the ghost problem were discussed. In the self-accelerating universe 
Htc = 1, there could be the enhanced symmetry that can eliminate the helicity-0 ghost if there were no 
spin-0 perturbation. Ref. [27] performed a gauge transformation 

A^^u = B^, + (V^V. + H^7i,.)X, X = (93) 

Then the action becomes 

Seff = ^j d'a^V^I - B^'-X^^iB) - H^B'^'^B^, + H^B^ 

- <^H{V^V. - 7^.n4 - 3if%^,)B''^ 

They propose to treat the spin-0 perturbation Lp as a, Lagrangian multiplier and perform quantization a 
la Nakanishi-Lautrup way in the QED. Effectively, this removes (/? from the spectrum and the enhanced 
symmetry ensures that there is no helicity-0 mode and no ghost (see Ref. [52] for a criticism on this proce- 
dure). However, this method cannot be applied if there is matter fluctuations with non-zero trace of energy- 
momentum tensor. This is because the trace of energy-momentum tensor breaks the enhanced symmetry 
and excite ^p as we saw in section IV. A. 

The other way to remove the ghost is to take into account the non-normalizable modes in the spectrum 
[27]. Then it is possible to eliminate the double pole from the amplitude that was the origin of the ghost. 
Instead, the amplitude has a single pole corresponding to the non-normalizable massless mode. However, it 
is found that this massless mode has an opposite sign for the amplitude compared with continuous massive 
states. This indicates that this massless mode becomes ghost-like. A very similar phenomena was found in 
the analysis of the shock wave analysis where the non-normalizable modes contribute a repulsive potential 
[26, 53]. 

Modifying the model 

Finally, in order to remove the ghost from the theory, we may have to modify the starting action. There 
are several attempts to extend the model [28, 54-57]. As we explained in section III.D, it is impossible to 
remove the ghost in two-brane model [28] . It was also shown that the introduction of Gauss-Bonnet term in 



19 



the bulk does not help [54] . Recently, a new model was proposed where the bulk is the solution in the normal 
branch but the junction condition is the one in the self-accelerating branch [56]. In order to achieve this, the 
5D Einstein -Hilbert action with a opposite sign to the conventional one is assumed to be localized near the 
brane. It remains to be seen whether this model can evade the ghost or not by studying the perturbations. 
The other approach would be consider higher co-dimensional branes [58] or intersecting branes [59]. It still 
remains to be seen whether there exists the self-accelerating universe and there is no ghost in these models. 
It is also important to study whether it is possible to embed the DGP model in string theory [60]. A UV 
completion of the model is necessary to study the fate of the ghost instability and it can guide us to make 
some simple modification of the model that will ameliorate the problems discussed in this review [48]. 
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